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HIGHLIGHTS 


• A scenario optimization model addressing weather uncertainty in BSC is proposed. 

• The modeling objective is to minimize the cost of biomass supply to biorefinery. 

• Field harvest work hours influence major cost-related decisions in the BSC system. 

• Yield of biomass is a crucial factor in determining the feasibility of BSC system. 

. Biomass storage method selected is dependent on the cost and dry matter loss during storage. 
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The aim of this study was to develop a scenario optimization model to address weather uncertainty in the 
Biomass Supply Chain (BSC). The modeling objective was to minimize the cost of biomass supply to 
biorefineries over a one-year planning period using monthly time intervals under different weather 
scenarios. The model is capable of making strategic, tactical and operational decisions related to BSC 
system. The performance of the model was demonstrated through a case study developed for Abengoa 
biorefinery in Kansas. Sensitivity analysis was done to demonstrate the effect of input uncertainty in 
yield, land rent and storage dry matter loss on the model outputs. The model results show that available 
harvest work hours influence major cost-related decisions in the BSC. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Presently, biofuels account for around 3% of the transport fuel of 
the world and is estimated to increase by 30% or more by 2050 (BP, 
2013). It has also been estimated that 30% of the petroleum-based 
energy consumption in the U.S. can be replaced with biomass, 
which will require roughly 1 billion tons of biomass (Perlack 
et al„ 2005). But, producing and supplying large quantities of low 
density biomass to the biorefineries is challenging. Feasibility of 
alternative energy sources such as biodiesel and hydrogen is also 
investigated but still remains largely in research phase. Major 
barriers preventing commercialization of cellulosic biorefineries 
is its production technology, limited biomass availability, and com¬ 
plex Biomass Supply Chain (BSC) system (Sharma et al„ 2013). It is 
estimated that biomass supply accounts for 20-30% of the ethanol 
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production cost, of which 90% is associated with logistical 
processes (Eksioglu et al„ 2009). 

Recently, there has been extensive focus by researchers on 
developing mathematical models for BSC design and management. 
Majority of the research works have developed mixed-integer lin¬ 
ear programming models with decision making capability ranging 
from strategic to operational-level (Sharma et al„ 2013). Com¬ 
monly used quantitative performance measure for BSC models 
are cost minimization or profit maximization (Gunnarsson et al., 
2004; Huang et al„ 2010; Tembo et al„ 2003; Zamboni et al„ 
2009). New approaches for modeling BSC are also introduced such 
as state-task-network (Dunnett et al„ 2007), spatially explicit 
(Zamboni et al., 2009), multi-stage (Huang et al., 2009), multi-ech¬ 
elon (Dal-Mas et al., 2011; Giarola et al., 2012), time-staged multi- 
commodity (An et al., 2011), multi-objective/multi-period (You 
et al., 2011; You and Wang, 2011), two-stage linear programming 
(Cundiff et al„ 1997), techno-economic system model (Svanberg 
et al., 2013) and two-stage stochastic programming (Chen and 
Fan, 2012) models. These approaches have increased BSC decision 
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making capabilities, and also address some of the critical issues 
and complexity associated with the BSC system. Studies by You 
and Wang (2011) and You et al. (2011) provide comprehensive 
assessment of the BSC with a focus on economic, environmental, 
and social impact of biofuel production. Majority of the models 
developed for BSC emphasize on considering sources of variability 
due to process and environment into the models for better man¬ 
agement (An et al., 2011; Huang et al., 2010; You and Wang, 
2011). Some of the models developed for BSC consider uncertainty 
in demand and price by formulating different scenarios (Akgul 
et al., 2010; Dal-Mas et al., 2011; Zamboni et al., 2009). The impact 
of weather uncertainty on BSC is crucial as it limits the amount of 
biomass supplied to a biorefinery and this variability is not consid¬ 
ered by majority of the models. 

Cundiff et al. (1997) developed a two-stage linear program¬ 
ming model for herbaceous biomass supply from 20 different 
farm locations to a centrally located biorefinery. The model 
determines monthly material flow and storage capacity expan¬ 
sion for each producer for four weather scenarios. However, 
the model does not capture complex BSC structure, and does 
not estimate number of equipment units required and storage 
treatments used. Therefore, to explicitly account for weather 
uncertainty and provide enhanced decision making capabilities 
for BSC, a scenario optimization model was developed in this 
present study. The component of weather uncertainty was 
incorporated into the model by estimating work hours available 
for harvesting biomass. The model provides decision about 
acres leased, material flow, number of harvesting units, in-field 
transportation units and transportation units purchased and 
rented, allocation of machinery units, and storage treatments 


used. The model also considers the technical and operational 
characteristics of the machinery units, and before-frost and 
after-frost harvesting of biomass. The specific objectives of the 

• To formulate a scenario optimization model for biomass supply 
to a biorefinery under weather uncertainty. 

• To develop a case study for switchgrass supply chain to the 
Abengoa Biorefinery (AB) at Hugoton, Kansas. 

2. Methodology 

2.1. Mathematical model 

A scenario optimization model was developed to minimize cost 
of biomass supply to a biorefinery considering harvest, transporta¬ 
tion, and storage costs. One year planning period with monthly 
time increments was considered. The model utilizes the yearly 
weather data to make the harvesting work-hour decision and each 
year was considered as a weather scenario. The daily weather data 
was obtained from the Oklahoma Mesonet for determining the 
work-hours available for harvesting (Mesonet, 2012). Each weath¬ 
er scenario was assigned equal probability of occurrence as the 
weather pattern was considered to be random and unpredictable. 
The soil moisture content, rain, snow, and daylight hours deter¬ 
mine the number of harvest work-hours available in a time period. 
The network structure consists of biomass source sites, storage 
sites, and a biorefinery site. The harvest unit consists of a self-pro¬ 
pelled windrower, a rake with a tractor and a large square baler 
with a tractor. The in-field transportation unit comprised of a bale 



Fig. 1 . Schematic of various components considered in the model. 









































B. Sharma et al./Bioresource Technology 150 (2013) 163-171 


165 


handler and a bale stacker. The road transportation unit consists of 
a semi-truck with trailer. Fig. 1 presents the schematic of different 
components considered in the model. Notations for indices, param¬ 
eters, and decisions variables used for the model formulations are 
summarized in Tables 1 and 2. 

2.1.1. Objective function 

The cost minimization objective function for the BSC system is 
represented by Eq. (1). The cost components for this economic 
objective function were biomass procuring, transportation, storage, 
rented and purchased equipment cost presented in Eqs. (2)—(13). 
The value of excess biomass sold to outside sources in a weather 
scenario (where more than the required biomass is produced) 
was subtracted from the total cost. Constraints for the model are 
presented in Eqs. (14)—(34). 

Minimize Cost = C^ + C a i + C fchu + C vchu + C /c9tu + 

+ C/ctu + Ctctu + Cmtu + Ccbos + Ccstor - P s fl (1) 


Table 1 

Model indices and input parameters notations and description. 




W 

Parameters 

ACRELAND,, 

BFtw 

CACRE,, 

CAPHUAF 

CAPHUBF 

CAPSTOR 

CAPTU 

CAPTUL 

COUTDEM 

CPROC 

CRHUBF 

CRHUAF 

CRNPTUL 

DEN, 

DMLOSS s 

FCHUAF 

FCHUBF 

FCSTORs 

FCTU 

FCTUL 

INVENDj 

INVSTARTj 

LEN,w 

PFL 

PROP,, 


PWSw 

VCHUAF 

VCHUBF 

VCTUIBj 

VCTUL 

VCTUSBi 

VCTUSIij 

WDHt„ 

YADJ, 

YIELD, 


Description 

Set of biomass source sites indexed by i 
Set of inventory sites indexed by j 
Set of land categories indexed by 1 
Set of storage method indexed by s 
Set of time periods indexed by t 
Set of weather scenarios indexed by w 
Description 

Acres of land under land category 1 in source site i (acres) 
Before-frost time period t in weather scenario w (0 or 1) 

Cost of leasing acres for biomass production under land 
category 1 in source county i ($/acre) 

Capacity of harvest unit after-frost (acres/hr.) 

Capacity of harvest unit before-frost (acres/hr.) 

Capacity of storage unit (tons) 

Capacity of road transportation unit (tons/load) 

Capacity of in-field transportation unit (tons/hr.) 

Cost of procuring biomass ($/ton) 

Rental rate for the harvest unit before-frost ($/day) 

Rental rate for the harvest unit after-frost ($/day) 

Rental rate for the in-field transportation units ($/day) 
Demand of biomass in a time period t (tons) 

Dry matter loss for storage method s (fraction) 

Annual fixed cost of harvest unit after-frost ($/year) 

Annual fixed cost of harvest unit before-frost ($/year) 

Fixed cost of storing biomass with storage method s ($/unit) 
Annual fixed cost of transportation unit ($/year) 

Annual fixed cost of in-field transportation unit ($/year) 

Road round trip hours (1/hr.) 

Ending inventory at inventory site j (tons) 

Beginning inventory at inventory site j (tons) 

Length of time period t in weather scenario w (hr./tp) 

Value of biomass sold to outside source/feedlot ($/ton) 
Proportion of land category 1 at source site i suitable for 
biomass cultivation (fraction) 

Probability of weather scenario w (fraction) 

Variable cost of harvest unit after-frost ($/acre) 

Variable cost of harvest unit before-frost ($/acre) 

Variable cost of transporting biomass from inventory site j to 
biorefinery site ($/ton) 

Variable cost of in-field transportation unit ($/ton) 

Variable cost of transporting biomass from source site i to 
biorefinery site ($/ton) 

Variable cost of transporting biomass from source site i to 
inventory site j ($/ton) 

Number of harvesting workdays available in a time period t 
and weather scenario w 

Yield adjustment factor for time period t (fraction) 

Yield of biomass for land category 1 (tons/acre) 


Biomass procuring cost (Cb P ) 

= (ttitw* w X CPROC x BF [W X PWS vv j 

+ (EEEEj^ 1 '* x cproc x o - Bftw ) x ( 2 ) 

Cost of acres leased (C a i) 

C a , = ^EEE acrelease,i x CACRE, t x PWS w j (3) 

Fixed cost (Cf c hu) and variable cost (C VC hu) of harvest units 

C fchu = (hubf x FCHUBF) + (ahuaf x FCHUAF) (4) 

C K hu = ^EEEE sab -fe» x VCHUBF x x PWS w j 

+ (eEEE^- x VCHUAF x (1 - BP W ) x PWS w j 

(5) 

Fixed cost (C/q/tu) and variable cost (C vc ift u ) of in-field transportation 
C fciftu = ntul x FCTUL (6) 

CmU = (EEEE W*» x VCTUL x BF <” x PWS ») 

+ (EEEE ta f»™ X VCTUL x (1 - Bftw) X PWS w j 

(7) 

Fixed cost ( C/ cn J and variable cost (C vctu ) of road transportation units 
C fctu = ntux FCTU (8) 

Cvcrn = (EEEE^W X VCTUSIij x PWS w ) 

+ ^EEE tssbitw X VCTUSBi X PWsSj 

+ (EEE^V x VCTUIBj x PWsSj (9) 


Variable cost (C vcru ) of rented harvest and in-field transportation units 
Cvcnx = ^EEE rn/!ub /™ x CRHUBF x BF m x WDAySJ 

+ (EEE rnhua -^ x CRHUAF x (1 - Bftw) x WDAYmJ 
+ (EEE renptu, *w x CRNPTUL x WDAYSJ 

(10) 

Cost of purchasing biomass from outside source (C C bos) 

C ch0 s = EE 0Ufdem tw X COUTDEM X PWS w (11) 

Cost of storing biomass with a specific storage method (C cstor ) 
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Ccstor = J2nstor s x FCSTOR s (12) 

Value of biomass sold to feedlot (P sfl ) 

p - lEEEE^ x PFL x pws w) (13) 


2.3.2. Constraints 

2.1.2.1. Supply constraints. Eqs. (14) and (15) presents constraints 
that ensure that the total acres harvested must not exceed the total 
acres leased for biomass production. The acres leased were equiv¬ 
alent to proportion of the total available acres in the biomass 
source site. The proportion of a land type available for biomass pro¬ 
duction in a biomass source site was estimated based on previous 
literature and can be varied depending on geographical area and 
land type. 

J2(sabf itw X BFt„ + saaf litw x (1 BFm)) 

= acrelease,, x PROP,, VI, i, w (14) 

acrelease,, V ACRELAND„ x PROP u VI, i (15) 

Eqs. (16) and (17) ensure that the total tonnage procured by the 
biorefinery equals the acres of biomass harvested multiplied by the 
yield and yield adjustment factor. The yield adjustment factor var¬ 
ies between 0 and 1. and adjusts the change in yield of biomass 
with progressing harvest season. 

tbfi itw = sabfitw x YADJ , x YIELD, x BF lw VI, i, t,w (16) 

tcifmw = saaf, itw x YADJf x YIELD, x (1 — BF m ) Vl,i,t,w (17) 


2.3.2.2. Balance constraints. The balance constraints are imposed on 
the biomass source site, biorefinery site, and inventory site. It was 
assumed that biomass was not stored at the biomass harvesting 
site. Therefore, the amount of biomass harvested at the biomass 
source site was equal to the total amount of biomass transported 
to the inventory site and the biorefinery site and the excess bio¬ 
mass if any sold to the feedlot (Eq. (18)). 

J2(tbf mw x BFtw + taf, iM x (1 - BFtw)) 

= tssb itw + ytssi i]M + ytsfl, itw Vi,t,w (18) 

At the biorefinery site, the balance constraint represents that 
biomass transported to the biorefinery site from the source site 
and inventory site plus the unmet demand of biomass must be 
equal to the demand of biorefinery (Eq. (19)). If there is unmet de¬ 
mand, biomass will be purchased from the outside sources. 

J i 

ytsibjM + ytssb llw + outdem [w = DEN t Vt, w (19) 

At the inventory site, the total biomass stored at inventory site 
in a time period is equal biomass transported from the source site 
to the inventory site in a time period plus biomass stored at the 
inventory site at the end of previous time period less any storage 
losses (Eq. (22)). The constraints (Eqs. (20) and (21)) allocate bio¬ 
mass transported to the inventory site from biomass source sites 
into different storage methods. Eq. (23) shows that the safety stock 
of biomass should be maintained at the inventory site at the end of 


planning horizon. The amount of biomass stored at inventory site 
was constrained by the storage unit capacity and number of stor¬ 
age units endogenously determined by the model (Eq. (24)). 


ytssiijtw = ytasip™ VJ,t,w (20) 

tsib jnv = yttiSj Stw Vj,t,w (21) 

stijstw = INVSTARTj + tasi jstv/ - ttiSj Stw + (if(t > 1) 

(sti M ^ )w x (1 - DMLOSS s , 0))) Vj, s, t, w (22) 

s% i,t i2,w >- INVENDj + youtinv M Vj, w (23) 

y S ti pM <= ns tor s x CAPSTOR Vs,t,w (24) 


2.1.2.3. Harvest units constraints. The total number of harvest units 
after-frost is equal to the number of harvest units used during the 
period before-frost along with additional units required to handle 
the heavier demands of added acreage for after-frost harvesting 
(Eq. (25)). 

huaf = hubf + ahuaf (25) 

Eq. (26) represent constraints that ensures the biomass acres 
harvested may not exceed the combined capacity of harvest units 
allocated to the biomass source site. The work hours available for 
harvesting biomass in a weather scenario and harvesting unit 
capacity, constraints the amount of biomass harvested at the 
biomass source site. Similarly, work hours available in a weather 
scenario and harvesting unit capacity constraints the amount of 
biomass transported during in-field transportation (Eq. (27)). Eqs. 
(28)-(33) represent constraints that ensure that the sum of the 
units (harvest and in-field transportation) used for biomass 
harvesting and handling at the biomass source sites may not 
exceed the total number of units determined by the model. 


(sabf, iav x BFm + saaf, itw x (1 - BFtw)) < 

= ( nhubf inv x CAPHUBF x BF lw + nhuaf ilw x CAPHUAF x (1 


-BF m ))xWDHm Vi, t, w 


(26) 

tbfiM x BF m + ta/iitw x (1 - BFtw 

) 


< nptul, x CAPTUL x WDH, w 

VI, i, t,w 

(27) 

nhubfi ,w = nphubf cw + mhubf itw 

Vi, t, w 

(28) 

nhuafjM = nphuafiM + mhuaf itw 

Vi, t, w 

(29) 

nptul icw = npurtul j[w + renptul itw 

Vi, t, w 

(30) 

ynphubf iM = hubf x BFm 'Vt, w 

(31) 

ynphuaf ,w = ahuaf x (1 - BF, V 

„) Vt,w 

(32) 

ynpurtul,M = ntul Vt, w 


(33) 
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variables 


[ description. 


Variables Description 

acreleaseu Total acres leased for biomass production under land category 1 in source site s 

ahuaf Number of additional harvest units purchased after-frost 

huaf Total number of harvest units after-frost 

hubf Number of harvest units purchased before-frost 

nhuafitw Total harvest units after-frost allocated to source site i in time period t and weather scenario w 

nhubfjtw Total harvest units before-frost allocated to source site i in a time period t and weather scenario w 

nphuaf itw Purchased harvest units after-frost allocated to source site i in a time period t and weather scenario w 

nphubfinv Purchased harvest units before-frost allocated to source site i in a time period t and weather scenario w 

nptulitw Total in-field transportation units allocated to source site i in a time period t and weather scenario w 

npurtul ltw Purchased in-field transportation units allocated to source site i in a time period t and weather scenario w 

nstor s Number of storage units of storage method types 

ntu Number of transport units purchased 

ntul Number of in-field transportation units purchased 

outdem tw Demand met from outside supply in time period t and weather scenario w 

outinvtw Inventory safety stock met from outside supply in time period t and weather scenario w 

renptulitw Rented in-field transportation units allocated to source site i in a time period t and weather scenario w 

mhuaf itw Rented harvest units after-frost allocated to source site i in a time period t and weather scenario w 

rnhubf itw Rented harvest units before-frost allocated to source site i in a time period t and weather scenario w 

saaf utw Acres of biomass harvested from land category 1 in source site i after-frost in time period t under weather scenario w 

sabfutw Acres of biomass harvested from land category 1 in source site i before-frost in time period t and weather scenario w 

sti jstw Tons of biomass stored at inventory site j under storage method s in time period t and weather scenario w 

tafutw Tons of biomass from source county i and land category 1 after-frost in time period t under weather scenario w 

tasi jstw Tons of biomass transported from source site i to inventory site j allocated to storage method s in time period t and weather scenario w 

tbf utw Tons of biomass from source county i and land category 1 before-frost in time period t under weather scenario w 

tsfliit„ Tons of biomass biomass harvested from land category 1 in source site i sold to feedlot in a time period t and weather scenario w 

tsibjt„ Tons of biomass transported from inventory site j to biorefinery site in time period t and weather scenario w 

tssb itw Tons of biomass transported from source site i to biorefinery site in time period t and weather scenario w 

tssiijtw Tons of biomass transported from source site i to inventory site j in time period t and weather scenario w 

ttiSjstw Tons of biomass transported from inventory site j to biorefinery site allocated to storage method s in time period t and weather scenario w 


2.1.2.4. Road transportation units constraints. This constraint en¬ 
sures that the biomass tonnage transported from the biomass 
source site to inventory site, inventory site to biorefinery site and 
source site to biorefinery site may not exceed the total number 
of transportation units (Eq. (34)). The transportation unit capacity 
is constrained by the length of time period and number of trips 
made by transportation unit in a time period. 

£ J i J 

J2J2 tss im« + J2 tssb itw + w 

< CAPTU X LEN [W X htrip Vt,w (34) 


2.2. Case study description 

The case study was developed for the Abengoa Biorefinery (AB) 
at Hugoton, Kansas (KS). AB considers switchgrass to be a potential 
future biofeedstock and therefore, it is considered as the biomass 
used by biorefinery for this case study analysis. In addition, switch- 
grass has a wide harvest window of eight months starting from 
July to February, which results in only four months of non-harvest¬ 
ing season. Twelve weather scenarios from 1998-2010 with equal 
probability of occurrence of 0.0831(1-^12) were considered. The 
50 mile radius was defined as the feedstock procurement area 
which includes Morton, Stanton, Grant, Haskell, Seward, and Ste¬ 
vens counties in KS and Texas County in Oklahoma (OK) (DOE/ 
EIS, 2010). For the operation of the biorefinery a daily demand of 
800 dry tons of switchgrass was considered with a nominal pro¬ 
duction schedule of 350 days per year (DOE/EIS, 2010). Since 
switchgrass is not commercially available in the feedstock procure¬ 
ment area, an assumption was made that the farmers would con¬ 
vert their land to switchgrass production. It was assumed that 
non-irrigated harvested cropland, permanent pasture, and expiring 
and not re-enrolled CRP acres 2011 and 2012 will be used for 
switchgrass production. The proportion of cropland and pasture- 
land used for switchgrass production in the source sites was 


adapted from the Billion-Ton Update report (Perlack and Stokes, 
2011 ) and was 25% and 50% for non-irrigated cropland and pas- 
tureland, respectively. It was also assumed that up to 50% of CRP 
land can be used for switchgrass production (Khanna et al„ 
2011 ). Yield of switchgrass decrease with maturity and the model 
assumes a yield loss factor of 5% per month from September 
through February (Epplin et al„ 2007). The AB anticipates an aver¬ 
age yield of 1.5-3.0 dry tons per acre for switchgrass from the soil 
types and climatic patterns present in the region (DOE/EIS, 2010). 
The yield of switchgrass was assumed to be two dry tons per acre 
for the base scenario. The farmers supplying biomass to AB were 
given three payment options (a fixed price option, a lower base 
price with revenue sharing and a lower base price with compensa¬ 
tion for nutrient replacement) for biomass cost per ton. The fixed 
price option of $15 per dry ton was used for this study (Bevill, 
2011 ). The cost of biomass purchased from outside source was as¬ 
sumed to be $90 per dry ton, which was 50% higher than the $60 
per dry ton delivered cost of biomass (Perlack and Stokes, 2011) 
and cost of selling excess biomass produced in a weather scenario 
to feedlot was assumed to be 25% of dollar per ton cost of biomass. 
Five inventory sites were selected based on road accessibility and 
the distance from the biorefinery site of about 30 miles. ArcGIS 
10 software (ESR1. ArcGIS 10.1) was used for selection of the inven¬ 
tory sites. One year safety stock (level of extra stock of biomass 
maintained at inventory sites to minimize the risk of stakeouts 
or not having biomass available to meet the biorefinery demand) 
of biomass at the inventory sites was maintained. The ending 
inventory was assumed to be equal to beginning inventory. The 
distances between the sites was calculated using the great circle 
distance formula. 

Two storage methods were considered: unprotected bales on 
the ground and tarped bales on gravel with dry matter loss of 
0.33% for each method. The dry matter loss values for each storage 
method were obtained from experiments conducted in the 
Oklahoma Panhandle for storing square switchgrass bales 
(E. Miller, 2012, Dry matter loss for square bales stored in 
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Panhandle, Oklahoma. Personal communication, Oklahoma State 
University). The reason for same dry matter loss values was 
attributed to extremely dry weather conditions during the storage 
period (2010-2011). The size of a storage bale stack was limited to 
100 tons and each stack was considered to be separated by at least 
20 ft., as required by the International Fire Code (ICC, 2003). 
Switchgrass can be harvested before-frost and after-frost. The first 
freeze in Oklahoma usually occurs in late October or early 
November (Koss et al., 1998; Sharma et al„ 2011). A Matlab 
program (MathWorks®. Ver., R2009b.) was developed to determine 
the time periods before-frost and after-frost. A day was considered 
to be before-frost if the minimum temperature was above 32°F. If 
60% of the days in a month met the criteria, the month was consid¬ 
ered to be “before-frost” for classification purposes. 

2.3. Machine performance 

2.3.1. Harvest unit (windrower, rake and baler) 

Biomass such as switchgrass can be mowed using standard hay 
equipment (Christensen and Koppenjan, 2010). A 16-ft. self-pro¬ 
pelled windrower was assumed to be used for harvesting switch- 
grass. The working throughput capacity of a self-propelled 
windrower, rake and large square baler is approximately 48, 39, 
and 21 dry tons per hour, respectively (M. Flerron, 2012. Through¬ 
put of harvest equipment. AGCO. Personal communication). The 
windrower operator adjusts speed of equipment to achieve the 
working throughput capacity under different yield conditions. 
The operator decreases the speed of equipment if the yield of 
switchgrass is high and vice versa. The speed adjustment was 
made by using Eq. (35) (Sokhansanj et al., 2008). If the adjusted 
speed of the windrower was outside the windrowers typical speed 
range provided by ASABE EP497.5 standards (ASABE, 2006), then 
the upper limit of the speed range was considered for the 
calculations. 

S = S b (y) (35) 

where, S = Adjusted speed (mph), S b = Optimum speed for the opti¬ 
mum yield (mph), Y b = Optimum yield (tons/acre), Y = Yield (tons/ 

A 13 ft. rake powered by a 75 hp tractor was considered for the 
present study. Under low yield conditions varying between one to 
two dry tons per acre, the rake gathers switchgrass from two sep¬ 
arate windrows into one windrow. The biomass was assumed to be 
packaged into 3ft. x 4ft. x 8ft. large square bales which allow bet¬ 
ter fit into the semi-trailer truck and maximize load per trip. 

Kemmerer and Liu (2012) found that the harvest system effi¬ 
ciency can be increased by balancing the field capacities of all 


equipment. The baler had the lowest material capacity in compar¬ 
ison to other operations such as raking, accumulating the bales in 
the field, loading the truck and stacking the bales in storage. The 
harvest unit was assumed to be integrated set of equipment con¬ 
sisting of a self-propelled windrower, a rake with a tractor and a 
large square baler with a tractor. The equipment in the harvest unit 
was matched to provide an Effective Field Capacity (EFC) of the 
large square baler. A harvest unit with an EFC of 25, 21, and 12 
acres per hours at the yield of one, two, and three dry tons per acre, 
respectively was considered (Table 3). 

2.3.2. In-field transportation and road transportation units 

The bale stacker collects, transports, and stacks bales at the side 
of the field. The Stinger Stacker 6500 handles 12 bales (3 ft. x 4 ft. 
or 3 ft. x 3 ft.) per load. The bale handler is used to load and unload 
bales onto semi-trailer trucks. The EFC of bale stacker and bale 
handler was assumed to be 34 and 32 tons per hour, respectively. 
It was assumed that a 52 ft. semi-truck trailer (21 tons of biomass 
per load) would be used for transportation of the biomass to the 
biorefinery. 

2.4. Economics 

The cost of machinery was broken down into two categories, 
the fixed cost and the variable cost. The proposed model considers 
procuring, harvesting, transportation, and storage costs. The har¬ 
vesting and transportation costs included the cost of harvest units, 
in-field transportation units, and road transportation units. The 
fixed cost of equipment was calculated using capital recovery 
method which included depreciation and considers wear, obsolesc- 
es, and age of machine. The variable cost was dependent on the use 
of machinery and comprised of repairs and maintenance (R&M), 
fuel and lubrication, and labor cost. 

2.5. Available harvesting work hours 

Switchgrass can be harvested for eight months starting in July 
until February of next year (Mapemba, 2005). The weather condi¬ 
tions determine days available for agricultural operations such as 
planting, cultivation, harvesting, and baling. Weather not only 
affects the time available for the field operation, but also the effi¬ 
ciency with which the operation is performed (Khanchi et al., 
2013). Estimation of work-hours is crucial for making decisions 
for machinery management and farm planning (Rounsevell and 
Jones, 1993). The methodology used for estimating harvest 
work-hours considered weather parameters such as rainfall, snow, 
temperature, and wind speed. The Oklahoma Mesonet daily weath¬ 
er data for Steven County, KS and the adjacent counties (Morton 


Table 3 

Effective field capacity of equipment in a harvest unit. 
Equipment Working width (ft.) 

Yield: 1.00 dry tons/acre 
Self-propelled windrower 16 

Rake 13 

Large square baler* 32 

Yield: 2.00 dry tons/acre 
Self-propelled windrower 16 

Rake 13 

Large square baler* 32 

Yield: 3.00 dry tons/acre 
Self-propelled windrower 16 

Rake 13 

Large square baler 16 


Speed (mph) 


4.5 


Efficiency (%) EFC (acres/hour) Number of equipment in a harvest unit 


At low yield of switchgrass such as 1.00 and 2.00 dry tons per acre, baler working width is 32 ft. At low yields, two windrows will be brought together to get enough 
material for baler pick-up. 
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County, KS; Grant County, KS; Stanton County, KS; Seward County, 
KS, and Texas County, Oklahoma (OK) was used for analysis. The 
Oklahoma Mesonet provided Fractional water Index (FWI) values 
which correspond to the soil moisture content. FWI is a unitless va¬ 
lue, ranges from zero (very dry soil) to one (very wet soil). The cri¬ 
terion used by researchers to determine the number of workdays is 
based on comparing soil moisture content to a percentage of field 
capacity depending on the type of soil (generally varies from 90- 
99.5%). In the present study, the criterion used for determining a 
non-workday was based on FWI values greater than 0.8. The rain¬ 
fall criterion used for estimating the number of workdays lost was: 
0,1, 2, 3, and 4 workdays lost for rainfall range of 0.00-0.05, 0.06- 
0.19, 0.20-0.49, 0.05-0.99 and >1.00 inches, respectively (Enz 
et al., 1995). The work hours available per day were calculated 
using Daylight Hours Explorer software (DaylightHoursExplor- 
er008. 2010). The available harvest hours were assumed to be on 
80 percent of the daylight hours. Twenty percent of the harvest 
work hours were assumed to be lost due to the machinery break¬ 
down, unavailability of labor, and unforeseen events. 

3. Results and discussion 

The scenario optimization model was developed to minimize 
cost of the biomass supply to the biorefinery. The model deter¬ 
mines the system cost, material flow, and number of machinery 
units purchased or rented. The optimization model was solved 
using Xpress-IVE version 1.22.04 (Fair Isaac Corporation ®2001- 
2012). The model was first analyzed for the base scenario (using 
parameters values described in Section 2.2-2.4) and later sensitiv¬ 
ity analysis was done on yield, storage dry matter loss, and land 
rent. 

3.1. Harvest work hours 

The weather patterns during each weather scenario determined 
the number of suitable work hours for harvesting switchgrass in a 
certain time period. The actual weather conditions cannot be pre¬ 
dicted in advance: therefore, the decisions about the purchase of 
the machinery were made based on the past weather data. This 
might result in excess machinery capacity during some years and 
insufficient machinery capacity during other years. In the present 
study, twelve years of weather data was considered beginning 
from 1997 to 2010 with each year considered as different weather 
scenario. Fig. 2 represents the work hours available for harvesting 
switchgrass in each weather scenario and time period. The average 
monthly harvest work hours available in each weather scenario 
were found to be statistically different (P = 0.0004). It was ob¬ 
served that the average work hours available for harvesting 
switchgrass for all harvest time periods varied from a minimum 


of 185 h to the maximum of 237 h in the 2009-2010 andl999- 
2000 weather scenario, respectively. The harvesting units were 
purchased/rented and deployed to the switchgrass source sites 
by the biorefinery. Therefore, it was assumed that the harvesting 
operation will be done on a rigid schedule utilizing maximum har¬ 
vest hours available to minimize the cost of the biomass supply to 
the biorefinery. 

3.2. Base scenario model results 

As estimated by the model, for the base scenario the number of 
harvest, transportation, in-field transportation and storage stack 
units required were 6, 23, 7 and 3149, respectively with total cost 
of biomass supply system as $240,141, 00. It was estimated that 
the investment to the biorefinery was highest (50%) for the storage 
stack units. The storage stack unit is defined as a stack consisting of 
180 bales. In the present model, it was assumed that a safety stock 
of a one-year switchgrass supply was stored at the inventory sites 
with bales covered with tarp and placed on gravel. The safety stock 
of switchgrass maintained at the inventory sites added a signifi¬ 
cant cost to the overall cost of biomass supply system. The factors 
considered for making decisions on purchasing equipment are har¬ 
vest work hours, yield of biomass, and cost of purchasing and rent¬ 
ing machinery. Twenty-three road transportation and five harvest 
units units account for 28% and 17%, respectively of the total cost of 
investment in purchased units by the biorefinery. In the base sce¬ 
nario, the yield of switchgrass was assumed to be two dry tons per 
acre. At low yields of biomass, rake operation is required through¬ 
out the harvesting season to gather enough switchgrass into a 
windrow for baler pick-up. Therefore, in the base scenario it was 
assumed that the harvest unit before-frost is equivalent to the har¬ 
vest unit after-frost. The in-field transportation units cost was 
found to be relatively low compared to the cost of other units. 

3.2.1. Acres leased and harvested for switchgrass production 

From the model output, it was observed that permanent pasture 
and expired/not re-enrolled CRP land acres were selected to be 
leased for switchgrass production in the biomass source sites. As 
yield of switchgrass was assumed to be two dry tons per acre for 
all land categories, the decision for land category acres contracted 
for switchgrass production was based on rental rate of land and the 
weather scenario pattern. The land rental values considered were 
$60, $40 and $55 per acre for non-irrigated cropland, permanent 
pastureland, and expired/not re-enrolled CRP land, respectively. 
Primarily, permanent pastureland was contracted for switchgrass 
production. In Seward, Stanton, Haskell, and Texas County, the per¬ 
manent pasture land contracted for switchgrass production was 
89%, 34%, 75% and 2%, respectively. Approximately, 5% of the total 
expired/not re-enrolled CRP land acres were contracted in Stevens 
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County by the biorefinery. The model assumes that demand of 
switchgrass by the biorefinery will always be satisfied. Therefore, 
the acres leased are sufficient to provide a regular and continuous 
supply of switchgrass based on yield, rental rate for land category, 
harvest, transportation and storage cost, and the weather scenario. 
All the contracted/leased acres were harvested in each weather 
scenario. The acres harvested in each time period of a weather sce¬ 
nario were dependent on harvest work hours available in the time 
period. The quantity of switchgrass available in each time period 
and weather scenario was consistent with the acres of switchgrass 
harvested. 

3.2.2. Allocation of total harvest units to switchgrass source counties 

The model allocates the harvest and in-field transportation 

units to the source sites depending on the quantity of biomass to 
be harvested at those source sites. Table 4 gives an example of 
the model output for a year 2009-2010, which shows allocation 
of harvest units during the harvest season to the switchgrass 
source counties. The model takes these decisions based on the 
work hours and the acres available for harvest from each site. Six 
harvest units were purchased by the biorefinery and additional 
capacity, if required for harvesting was achieved by renting the 
harvest units. It was observed that the weather scenario with least 
work hours rented the maximum harvest units. The 2009-2010 
weather scenario had a total of 1476 work hours (minimum work 
hours among all weather scenarios) and therefore, rented five addi¬ 
tional harvest units to meet the capacity. A similar trend was ob¬ 
served for the allocation of in-field transportation units to the 
switchgrass source sites. 

3.2.3. Transportation of switchgrass 

The model makes the decision for transportation of switchgrass 
based on two options: the feedstock can go to an inventory site be¬ 
fore going to the biorefinery or it can directly go to the biorefinery. 
Depending on the demand of switchgrass, the quantity can be re¬ 
moved from the inventory sites in any of the time increments. 
For the time period with less harvest work hours lower quantities 
of biomass were transported to the biorefinery site. The demand 
for these time periods was satisfied by removing switchgrass from 
inventory sites. Switchgrass was stored at inventory sites to meet 
the demand of 120 days of non-harvesting period from March to 
June of a weather scenario. It was observed that switchgrass was 
transported to inventory sites starting from the first time period 
to build the inventory for non-harvesting season. In some time 
periods, lower quantities or no switchgrass was transported to 
the inventory sites. This can be attributed to lower number of har¬ 
vest work hours in a time period. In case the quantity of switch- 
grass transported to inventory sites was not enough to meet the 
demand of non-harvesting time periods/months, switchgrass was 
purchased from outside source at a higher price. It was observed 
that in each time period and weather scenario, the total quantity 
of switchgrass harvested at source sites was equal to the sum of 
the quantity of switchgrass transported to the biorefinery site 
and inventory sites indicating no discrepancies in the model. 

3.2.4. Storage of switchgrass 

Switchgrass was stored and accumulated at the inventory sites 
to meet the demand of the biorefinery for non-harvest time peri¬ 
ods and maintain safety stock of switchgrass at inventory sites. A 
safety stock for one year supply of switchgrass (280,000 dry tons) 
was maintained at inventory sites for beginning (July) and ending 
(February) time periods of a weather scenario. Stored switchgrass 
was also used to meet the demand of the biorefinery for time per¬ 
iod with less harvest work hours, and also account for switchgrass 
losses during storage. In all weather scenarios, switchgrass was 
transported and accumulated at inventory sites from July-Febru- 


Allocation of total harvest units to switchgrass source counties in 2009-2010 weather 
scenario. 


Jul. Aug. Sep. Oct. Nov. Dec. Jan. Feb. 


Texas Co. (OK) 0 

Morton Co. (KS) 2 

Stanton Co. (KS) 0 

Stevens Co. (KS) 1 

Grant Co. (KS) 0 

Haskell Co. (KS) 0 

Seward Co. (KS) 3 


0 

2 

0 


0 


ary. From March -June time period, switchgrass was removed from 
inventory sites to meet the demand of the biorefinery for non-har¬ 
vesting periods. The ending inventory in the June time period was 
maintained at 280,000 dry tons. In weather scenario 2000-2001, 
switchgrass was removed from inventory sites in the December 
time period because less harvest work hours were available in 
the December and enough switchgrass was not harvested to meet 
the demand of biorefinery. The storage method selected was the 
uncovered bales placed on the ground. The model makes decisions 
on storage methods selected for storing biomass at the inventory 
sites based on dry matter loss and cost of the storage method. As 
mentioned earlier, the dry matter loss was same for both the stor¬ 
age alternatives and cost of storing bales unprotected on the 
ground was less compared to storing bales with tarp placed on 
gravel pad, therefore uncovered method of storage was selected. 


3.3. Sensitivity analysis 
3.3.1. Yield of switchgrass 

The profit for the biorefinery is influenced by several factors 
such as feedstock type, yield, harvest work hours, biorefinery size, 
biorefinery location, storage losses, and transportation cost. Yield is 
one of the most crucial factors considered for the selection of bio¬ 
mass for bioethanol production and substantially changes land 
usage, cost, and machinery requirement to meet demands of the 
biorefinery. For investigating the influence of yield on model re¬ 
sults, the scenarios of one, two, and three dry tons per acre were 
used as input to the model. It was observed that at the low yield 
of one dry tons per acre, all the switchgrass required by the biore¬ 
finery was purchased from outside source. In a real world situation, 
even at low yield of biomass, the biorefinery will harvest the acres 
contracted for the biomass production. The present model is devel¬ 
oped for a one year planning horizon. Therefore, under a low yield 
of biomass, the least cost option was to procure switchgrass from 
outside sources. Among the three yield scenarios, the cost of the 
switchgrass supply to the biorefinery was least for the three dry 
tons per acre yield scenario. The difference in cost of switchgrass 
supply for two and three dry ton per acre yield scenarios was 
$1,249,500. At three dry tons per acre yield, the harvest unit com¬ 
prised of one self-propelled windrower, two large square bale, one 
rake and two tractors with the Effective Field Capacity (EFC) of 
12 acres per hour. When the yield was two dry tons per acre, the 
harvest unit comprised of two self-propelled windrower, one large 
square bale, two rake and two tractors with the EFC of 20 acres per 
hour (Table 3). The difference in cost of harvest units at two and 
three dry tons per acre yield was $6,035. The reason for less varia¬ 
tion in biomass supply cost of two and three dry tons per acre yield 
scenario could be attributed to large difference in effective field 
capacity and low difference in cost of harvest units for the yield 
scenarios. Fewer acres were contracted at the yield of three dry 
tons per acre, as the demand was met due to higher yield of 
switchgrass. 
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3.3.2. Land rent and dry matter storage loss 

Land rent is a crucial factor in determining the willingness of 
farmers to allocate their land to biomass production. As land rent 
increases the farmers will be more willing to grow dedicated 
bioenergy crops. In the present study, the effect of increase in land 
rental values by 15% was evaluated. Two scenarios were consid¬ 
ered. Scenario-1 (base scenario) with the land rental values for 
non-irrigated cropland: $60 per acre, permanent pasture: $40 per 
acre, and expired- CRP: $55 per acre. Scenario-2 with the land 
rental values for non-irrigated cropland: $69 per acre, permanent 
pasture: $46 per acre, and expired- CRP: $63 per acre. It was ob¬ 
served that land rental rate had a significant impact on decisions 
for leasing acres for switchgrass production. It was observed that 
only permanent pasture land was contracted in scenario-2. In sce¬ 
nario-1 permanent pasture land and expired/not re-enrolled CRP 
land acres were contracted by the biorefinery. The other variables 
did not vary significantly. 

In the base scenario, the dry matter loss was the same for both 
storage alternatives and the cost of storing bales uncovered was 
low. Therefore, the model selected the uncovered storage method 
as the best alternative for storing switchgrass. $ensitivity analysis 
was also performed to evaluate the effect of dry matter loss on 
allocation of the method selected for switchgrass storage. For the 
present study, the base scenario (scenario-1) was compared to sce¬ 
nario-2 with dry matter loss for unprotected (bales uncovered and 
placed on the ground) and protected (bales covered with tarp and 
placed on gravel) storage methods as three (0.66%) and four 
(1.33%) times higher than the value considered in the base sce¬ 
nario, respectively. The results indicate that as dry matter loss in¬ 
creases, the acres contracted to meet the demand of biorefinery 
also increases. One extra harvest unit was also purchased for har¬ 
vesting extra acres in the case of high dry matter loss scenario. 

4. Conclusion 

A scenario optimization model considering weather uncertainty 
was developed. Harvest work hours influenced major cost-related 
decisions in the B5C. The method of biomass storage was depen¬ 
dent on the cost and dry matter loss during storage. Yield of bio¬ 
mass was found to be a crucial factor in determining the 
feasibility of B5C system. The model captures weather uncertainty 
and its influence on purchasing and deployment of assets. The 
direction for future research will be to consider different types of 
biomass feedstocks, conversion processes, and environmental im¬ 
pact of biofuel production into the model. 
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